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Recently, superconductivity in PuCoGas was discovered. It has the same crystal structure 
as CeMIns(M =Ir, Co, Rh), which are often refered to as Ce-115's. The electron correlation 
in PuCoGas is estimated to be weak compared with Ce-115's, and the filling number of 
electrons is considered to be far from 0.5/spin in the band which plays an important role 
in realizing the superconductivity. Nevertheless, the superconducting transition temperature 
T c in PuCoGas is almost by an order of magnitude higher than that in Ce-115's. In order 
to explain the superconductivity with high T c , we adopt the periodic Anderson model and 
calculate T c by solving the Dyson-Gor'kov equation derived by the third order perturbation 
theory with respect to U. By this calculation, we indicate that the superconducting state of 
PuCoGas is a <i-wave pairing state, and show that the good location of two Fermi surfaces 
results in the high T c in PuCoGas. 

KEYWORDS: unconventional superconductivity, periodic Anderson model, heavy Fermion, plu- 
tonium 



1. Introduction 

Recently, superconductivity in PuCoGas was discovered by Sarrao, et al. It has very 
high transition temperature (T c = 18. 5K). This value of T c is higher than that in any other 
isostructual superconductors, such as CeMIns (M =Ir, Co, Rh). Superconducting transition 
in CeCoIns and Celrlns occur at ambient pressure at T c = 2.3K and 0.4K, respectively. 2,3 But 
CeRhlns becomes superconducting only under pressure with T c = 2. IK. 4 Although NpCoGas 
and UMGa5 have the same HoCoGas-type crystal structure, superconductivity has never been 
reported in these materials. From now on, we refer to these HoCoGas-type compounds as '115'. 

First of all, let us consider Ce-115's. In the phase diagram of Ce-115's, antiferromag- 
netic(AF) state and superconducting state are adjacent to each other. 5 Moreover, the magnetic 
field dependence of thermal conductivity 6 and T 3 -behavior of nuclear spin relaxation rate in 
the superconducting state 7 show the existence of line-node gap. Band calculation shows that 
the Fermi surfaces of Ce-115's are quasi-two-dimensional. From these facts, Ce-115's have been 
considered to be unconventional quasi-two-dimensional d-w&ve superconductors induced by 
antiferromagnetic spin fluctuations (AFF). 8 Using Hubbard model, Nisikawa et al. explained 
that the superconductivity of Ce-115 has <i x 2_ y 2 symmetry. 9 
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Then, what is the mechanism of superconductivity in PuCoGas? Experimental facts such 
as the Curie-Weiss behavior in magnetic susceptibility at T > T c , T 135 -behavior in electric 
resistivity at T c < T < 50K, and power-law behavior in specific heat at T < T c are reported. 1 
Moreover, band calculations by Maehira et al. 10 and Opahle et al. 11 show that the Fermi 
surfaces of PuCoGas are quasi-two-dimensional just like Ce-115's. These facts imply that the 
superconductivity in PuCoGas is an unconventional d-wave superconductivity with magnetic 
origin, just like Ce-115's. Of course, there are differences between PuCoGas and Ce-115's. One 
is the value of T c , and another is the strength of electron correlation. Specific heat coefficients 
7 = C/T\ T=Tc , which are proportional to the renormalized electron mass m*, are 290, 400 
and 750 mJ/mol • K 2 for CeCoIns, CeRhlns and CeLTns, respectively. On the other hand, 
7 for PuCoGas is 77mJ/mol • K 2 . This means that the mass of electron is not so enhanced 
and the electron correlation is modest in PuCoGas compared with Ce-115's. The theoretical 
specific heat coefficient 7band estimated from the band calculation 10 is 16.9mJ/mol -K 2 . Thus 
we can estimate the mass enhancement factor m* /m = j/-fband i n PuCoGas at 4.5. This value 
is rather lower than that in Ce-115's, which is more than 10. 

From Ref. 10 we can see that in PuCoGas there exist no bands which are near the half- 
filled. From this fact, it seems difficult to explain T c in PuCoGas which is high almost by 
an order of magnitude compared with Ce-115's. Let us see the Fermi surfaces of PuCoGas in 
band calculations. 10,11 The 16th band and the 17th band in Ref.10 have Fermi surfaces. Here 
we ignore the Fermi surfaces of the 15th band and the 18th band since they are very small. In 
these situations, the following points are important. Since PuCoGas nas the two main Fermi 
surfaces, the effective density of state at the Fermi energy becomes large. Furthermore, if two 
Fermi surfaces are well located, the effective correlation for antiferromagnetic fluctuation is 
strengthened as shown below, even though 7 is still not so large. In this paper we point out 
that this leads to relatively high T c for PuCoGas. 

2. Periodic Anderson Model 

Let us introduce the following periodic Anderson model. 12 The Hamiltonian is 



4 




(i) 



(2) 



(3) 



Vfc = Vq — V\ cos k x cos k y . 



(4) 
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Here, t and t' denote the nearest and next nearest neighbor hopping terms of /-electron, 
respectively. t c and t' c denote those of conduction electron. Thus, e, and e°, are the dispersion 
of /-electron and conduction-electron, respectively. is the hybridization between /-electron 
and conduction-electron. We set these parameters so that the band structure and the Fermi 
surfaces of diagonalized bands reproduce those of the band calculation in Ref.10 and Ref.ll. 
Hereafter we use the following parameters: t c /t = 6.0, t'Jt = 1.8, Vo/t = 2.8, V\/t = 2.1, 
Hc/t = 0.8 and t'/t = 0.3. The total number of filled electrons ntat in the /-band and the 
conduction-band is 1.16 per spin (58% filled). The chemical potential hq a t temperature T is 
determined by 

^E(/K)+/(4))=— (5) 

k 

where / (x) = (e x + l)" 1 is the Fermi distribution function. The unperturbed term of the 
Hamiltonian of Eq.(l) is rewritten in the 2x2 matrix form as follows 




(6) 



(7) 

Here, the /-band and the conduction-band are hybridized by Vu , and then diagonalized into 
two bands: the band-1 and the band-2. Operators f^ a ^/^ ^ , [c^ ^ , (a)^ ^ and 
bfc a (b^ ^ are the annihilation(creation) operators of the /-band, the conduction band, the 
band-1 and the band-2, respectively. E\ and E 2 are the dispersions of the diagonalized bands. 
{Ei > E 2 ) 

We now can express the bare Green's function as follows 



Go(k) 



c 2 G! + s 2 G 2 
ac(Gi - G 2 ) 




(8) 
(9) 
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sc= , Vk = (10) 



f c 

' 2 - * + y k \ (11) 

s 2 = 1 - c 2 . (12) 



In Eq.(9), e n = (2n + 1) 7rT is a fermion Matsubara frequency. 
3. Calculation by TOPT 



In our numerical calculation, we divide the first Brillouin zone into 128 x 128 meshes and 
take 4096 Matsubara frequencies. To treat electron correlation effect, we need to approximate 
the self-energy terms. Among several ways of approximation, we adopt the third-order per- 
turbation theory (TOPT) with respect to U. Using TOPT, we can write the self-energy terms 
as 

Si(k) = ^J2 V A^k')G f (k'), (13) 

v 

V n (k, k>) = ^U 2 Xf0 (k - k!) + ^U 3 [ X % (k - k') + cj>% (k + k')] , (14) 



N /VJU v ' N 

(a) = - 

N 



k 



4>fo{q) = - r ^Y, G i^) G i^- k )- (I 6 ) 

k 

Here, we have introduced the abbreviation k = (k,e n ) and q = (q,uj n ). Note that uo n = 2nirT 
is a boson Matsubara frequency. The dressed Green's function G(k) in normal state is given 
by 

G(k) = G (k) + G (k)S (k) G(k), (17) 

where, 



c 



\ ^0 ^0 

" V G c f G c 

[si - <y M ) o 





(18) 
(19) 
(20) 
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(-71,-Jl) (^,-71) 

Fig. 1. The Fermi surfaces of the band- 1 (dashed curve) and the band-2 (solid curve). 



The shift of chemical potential 5p is determined by conservation of total electron number 

Y,(G* + G c -Gl-G%) =0. (21) 

k 

It is noted that the self-energy correction appears only in /-electrons, since the electron 
correlation U is taken into account only among /-electrons. 

Now, we can calculate E[ and E' 2 , which are the modified dispersions of E\ and E2 by 
including the self-energy correction, respectively. They are given by 



2 2 



e% + e f k + E?(k,e = 0)-6 fi ) ± J U+ £R(k,e = 0) - 6»- tfY + 4V* 



, (22) 



where S^(k,e) is the retarded normal self-energy, which is calculated by analytic continu- 
ation ie n — » £ from i7n (A;). The Fermi surfaces and the band structure of the diagonalized 
bands calculated from Eq.(22) are shown in Fig.l and Fig. 2, respectively. Let us introduce the 
scaling parameter to as one eighth of the width of the band-1 for the sake of easy comparison 
with usual calculations for a single band Hubbard model, which corresponds to Vu = in the 
periodic Anderson model. In figures except for Fig. 3, we rescale all energies (such as U and 
E[) by to (f° r example, shown as U/to and E[/to). Note that to is not equal to t, which we 
set to be unity. The ratio to/t is 3 ~ 3.5 and depends on U. When U/to = 4.1, to/t is about 
3.1. Fig. 3 shows the bare density of state without self-energy correction. In Fig. 3 energies are 
rescaled by t, not by to- The bare density of states of /-band in the periodic Anderson model 
pf (e) and that in single band model p° (e) are given by 

Pf(e) = -^IraJ2G f R (k,e), (23) 
k 

P° f (s) = --lmY,Gl R (k,e), (24) 
k 
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(tt,tc) (0,0) (ti,0) (71,71) (n,2n)(2n,2n) (71,71) 



Fig. 2. The dispersion of the band- 1 (upper) and the band-2 (lower). 



, respectively. Here, the unperturbed Green's functions Gq R and Gq R are calculated by analytic 
continuation ie n — > e from Gq (k) and [\e n — ejQ , respectively. In Fig. 3, the DOS in the 
periodic Anderson model becomes flat and there exist high and low energy tails in the DOS 
owing to the hybridization. These facts mean that the width of /-band is expanded compared 
to that in single band model. The expanded band width stabilize the Fermi liquid state. Thus 
our perturbation calculation can be valid owing to the expansion of the band width, even 
when U exceeds 8t, which is the band width in the single band model. 

We now calculate T C D In our model, the Coulomb repulsion works only between f- 
electrons, so we have only to account of /-electron Green's functions. Then, the anomalous 
self-energy S a (k) is given by 

S a (k) = S RPA (k) + i7 vcrt (k) , (25) 

Srpa (k) = ~ Yl i U + ( k + k ') + 2C/3 4o ( k + k ')] F ( k ') . ( 26 ) 

k' 

S vert (k) = -U 3 G o (X/O (k + hi) - cf> f0 (k + k{))G{ (k + h- k') F (k') 

^ ' k',ki 

~ 1/3 (jf) E G o (X/O (~k + h) - <P f o (~k + h))G f {-k + h- k') F (k') . (27) 

^ ' k',ki 

In superconducting state G$ (k) and the anomalous Green's function for /-electron F (k) 
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Fig. 3. The DOS for /-electrons in periodic Anderson model (solid curve) and in single band model 
(dashed curve). In periodic Anderson model, the band width is expanded and the DOS is flat 
compared to those in single band model. 



satisfy Dyson-Gor'kov equations. 13 

G f (k) = G{ (k) + G{ (k) El (k) G f (k) + G{ (k) S a (k) ^ (A;) , 

Ft (k) = G f (-k) Si (-k) Ft (k) + Cq (-k) K (-k) Gf (k) . 

The normal self-energy is calculated in Eq.(13). 
In the vicinity of T c , F (k) can be linearized as 

F(k) = \Gf (k)^ S a (k), 

G f (k) = G f (k) + G f (k) S( (k) G f (k) . 

Thus, S a (k) at T c is determined by the gap equation 



where, 



V a (k, k')=U + U 2 X0 (k + k') + 2U\l (k + k') 
+ 2f/ 3 ^Re^ G f (h) G f (k + k ± - k') [xo (k + h) -<f> Q {k + h)}. 



ki 



(28) 
(29) 



(30) 
(31) 

(32) 



(33) 



If we replace the left hand side of Eq.(32) by \E a (k), this equation can be considered 



7/12 



Full Paper 





as an eigenvalue equation with eigenvalue A and eigenvector S a (k). T c is the temperature at 
which the maximum eigenvalue reaches to unity. S a (k) represents the superconducting gap 
symmetry. Among several gap symmetries, the d x 2_ y 2 state possesses the maximum eigenvalue. 
Fig.4 shows the analytic continuation (k,e = 0) of the gap function S a (k). 

Now we consider the condition in which the perturbation theory in U is valid. We in- 
vestigate the behavior of retarded normal self-energy S^(k,e). At wavevector k = k\ = 
(97r/16, 97r/16), which is near the Fermi surfaces of the /-band, Re Zjf and ImA^ behave as 
shown in Fig. 5. From Fig. 5, we can see that near e = 0, Re = as and Im = f3e 2 . These 
facts show that the retarded normal self-energy behaves as the conventional Fermi-liquid. The 
perturbation calculation up to third order terms of U is confirmed to be valid. 
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(n,-n) 



Fig. 6. Quasi-particle mass enhancement factor z 1 (k) (U/t = 5.4, T/t = 0.007). The values of 
z^ 1 (k) on the Fermi surfaces of band-1 and band-2 are 4.5 ~ 5.0. 



Next, we investigate the mass enhancement factor for the /-band, z~ l (k) = 1 — 
dRe£n (k,e) jde. Fig.6 shows z^ 1 (k). If the negative contribution of the U 3 -term to z~ l (k) 
is large compared to the U 2 -term, then z~ l (k) becomes near or lower than unity. So the 
value of z^ 1 (k) in Fig.6 shows that the contribution of the C/ 3 -term to z^ 1 (k) is not so large 
and that perturbation calculation in U is valid. 

From Fig.l, the numbers of filled electrons in the band-1 and the band-2 are estimated at 
0.34/spin and 0.79/spin, respectively. In Fig.l, we can see that the shape of the Fermi surface 
of the band-1 is preferable to the d x 2_ y 2 superconductivity, just like high-T c cuprates. On 
the other hand, the low value of the electron number in the band-1 (0.34/spin) suppresses 
AFF, and lowers the peak of spin susceptibility near (tt, tt). The reduced AFF leads to the 
d-wave superconductivity with low T c . Surely the band-1 is important for superconductivity, 
but we cannot explain high T c of PuCoGas if we consider only the band-1. From these facts, 
we consider that the band-2 plays an important role in high T c superconductivity. 

Fig. 7 shows the spin susceptibility including the self energy correction. Line 1 is spin 

susceptibility Xfo (<?) = X2 (q) + X3 (<?) + Xi (<?) Ior the periodic Anderson model. Lines 2, 3 

and 4 are X 2 (q) = £ c 2 d (k) c 2 Gi (k + q), X 3 (?) = 2 Re £ c 2 d (k) s 2 G 2 (k + q) and X 4 (q) = 

k k 
£ s 2 C?2 (k) s 2 C?2 (k + q), respectively. There is a peak near (tt, tt) for \2 (q), which corresponds 

k 

to the spin susceptibility for the single band model composed of the band-1. Moreover, the 
exsistence of the coupled contribution between the band-1 and the band-2 (x3 (<?)) causes the 
higher peak of the spin susceptibility, which is stronger AFF in our model, compared with that 
in single band Hubbard model. The large susceptibility arises from nesting effects between 
two bands. The higher T c is originated from this effect. The calculated T c in the periodic 
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(0,0) (0,71) (71,71) q (71,0) 



Fig. 7. Spin susceptibility for the single band model and the periodic Anderson model (U/t = 
5.4,T/i = 0.007). 



Anderson model is shown in Fig 8. The lowest limit of temperature for reliable calculation 
is approximately T/to = 0.002. From Fig. 8 we can see that T c is relatively high even for 
small values of U, or weak correlation. We tried to calculate T c in single band models, but we 
could not get any finite value of T c in TOPT because the filling number is far from 0.5/spin. 
This indicates that T c in the single band model is very low compared with T c in the periodic 
Anderson model. We can see that in PuCoGas the d-wave superconductivity with high T c is 
realized even for the modest electron correlation. The modest electron correlation is consistent 
with the experimental facts. 1 

4. Conclusion 

In this paper, we have explained high T c of PuCoGas using periodic Anderson model 
and have shown that the superconductivity in this material is unconventional one with d x 2_ y 2 
symmetry. To obtain the results, the following point is important. By considering only the 
'main' band (the band-1), it seems impossible to explain high T c because of low filling. In 
this material, the existence of 'sub' band (the band-2) increases much the density of states 
at Fermi energy. Furthermore, nesting effects between the Fermi surface of the 'sub' band 
and that of 'main' band enhance AFF. These effects make T c higher even for relatively weak 
electron correlation. From the band calculation, W\ ~ 1 eV ~ 1.16 x 10 4 K. 10 Since to is defined 
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Fig. 8. T c of PuCoGas calculated by TOPT on the basis of periodic Anderson model. 



as one eighth of W±, to is approximately 1.5 x 10 3 K. Roughly estimated, the superconducting 
transition temperature T c = 18. 5K corresponds to about 0.012to- From this value of T c and 
Fig. 8, we can estimate the value of U /to at 4 ~ 4.5. 

This calculation was performed with the computer in Yukawa Institute of theoretical 
Physics. The authors are grateful to Dr.Y.Nisikawa for valuable discussions. 
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